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Abstract. The recently published E865 data on charged Kei decays and tttt phases are reanalyzed to 
extract values of the two S-wave scattering lengths, of the subthreshold parameters a and /3, of the low- 
energy constants £3 and I4, as well as of the main two-flavour order parameters: (uu) and F-^ in the limit 
rriu = md = 0 taken at the physical value of the strange quark mass. Our analysis is exclusively based on 
direct experimental information on tttt phases below 800 MeV and on the new solutions of the Roy equations 
by Ananthanarayan et al. The result is compared with the theoretical prediction relating 2ag — Soq and 
the scalar radius of the pion, which was obtained in two-loop Chiral Perturbation Theory. A discrepancy 
at the l-cr level is found and commented upon. 


1 Introduction 

Recently, a new measurement of decay and of the tttt 
phase shift difference 5° — dj has been published [1] by the 
Brookhaven F865 collaboration. New low-energy tttt scat¬ 
tering data constitute a rare event which has not happened 
since the last Geneva-Saclay experiment 25 years ago [2], 
and the corresponding determination of the isoscalar S- 
wave scattering length Og = 0.26 ± 0.05. The new exper¬ 
iment [1] improves the statistics by more than a factor 
of 10, and the outcome for Og is not only more accu¬ 
rate, but also points towards a smaller central value as 
expected by the standard version of Ghiral Perturbation 
Theory (yPT) [3]. It is crucial to know the two S-wave 
scattering lengths Ug and Og as precisely as possible and to 
avoid confusing their model-independent extraction from 
the data on one hand with theoretical yPT-based pre¬ 
dictions on the other hand. This is the main purpose 
of the present paper. The principal model-independent 
tools of our analysis are Roy equations [4] and their re¬ 
cent solution by Ananthanarayan, Golangelo, Gasser and 
Leutwyler (ACGL) [5]. In a suitable kinematical range, 
Roy equations represent a rigorous consequence of gen¬ 
eral properties of the scattering amplitude: analyticity, 
unitarity, crossing symmetry and asymptotic bounds. The 
AGGL solution uniquely fixes the three phases relevant at 
low energies, 6q{E),61{E),Sq{E) for E < Eq = 800 MeV 
in terms of the two scattering lengths Ug and Cg, and ex¬ 
perimental data above sg = Eq. In this way, sufficiently 
precise data on phase shifts for 2M,r < E < Eq can be 


converted into a model-independent determination of the 
two S-wave scattering lengths. Gonversely, the knowledge 
of ttg and Qq allows one to determine any other low-energy 
TTTT observables which might be needed for a theoretical in¬ 
terpretation of experimental results in terms of the chiral 
structure of QCD vacuum. Unfortunately, the combina¬ 
tion dg — 1^1 near threshold (which is measured in 
experiments) is not strongly sensitive to Og, for a given 
ttg, even if the former is allowed to vary over the whole 
Universal Band (UB). This can be seen explicitly, for in¬ 
stance, from the ACGL solutions of Roy equations. We 
will show (Sec. 3) that a model-independent and rela¬ 
tively accurate determination of both S-wave scattering 
lengths is nevertheless possible, if the existing Ke 4 data 
on (5g — dj are combined with the older production data 
by Hoogland et al. [6] and by Losty et al. [7], concerning 
the 1 = 2 S-wave for E < 800 MeV. Our result reads 
Og = 0.228 ± 0.012, Ug = —0.0382 ± 0.0038. Subsequently 
(Sec. 4), Roy equations will be used once more to convert 
the S-wave scattering lengths into the determination of 
the subthreshold parameters a and f3 (as well as Ai ... A4) 
introduced in Ref. [8] and under a different name (61 ... fog) 
in Ref. [9]. The expansion of the parameters a and /3 in 
powers of quark mass converges more rapidly than in the 
case of scattering lengths. The knowledge of subthresh¬ 
old parameters is therefore used to extract the low-energy 
constants (LEC’s) £3 and £4, as well as the main two- 
flavour order parameters (Secs. 5 and 6). All this anal¬ 
ysis is based on existing solutions of Roy equations as 
given by ACGL in Ref. [5]. In parallel, we construct an 
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extended solution of Roy equations, incorporating the de¬ 
pendence on the value of phase shifts at the matching 
point = 800 MeV. This allows the control of the propa¬ 
gation of errors arising from the intermediate energy data, 
which was not possible using the published ACGL solution 
of Ref. [5]. Finally, our model-independent determination 
of scattering lengths and other low-energy parameters is 
compared with the theoretical prediction of the tight cor¬ 
relation between 2 ag — bog and the scalar radius of the 
pion (r^)s [10]. If the latter prediction is combined with 
the E865 data alone [1], one obtains Og = 0.218 ± 0.013, 
Og = —0.0449 ± 0.0033 assuming the value of the scalar 
radius (r^)s = 0.61±0.04 fm^. Our simultaneous fit to the 
Kei [1,2] and to the high-energy 7 r+ 7 r+-production data 
[6,7] suggests a discrepancy at the 1-cr level with the two- 
loop yPT prediction for 2aQ — Sog [10]. We point out that 
the treatment of symmetry breaking 0(p^) corrections in 
the scalar channel might be at the origin of this discrep¬ 
ancy (Sec. 7). 


2 Theoretical motivation 

Before we proceed further it may be useful to summarize 
briefly the main theoretical motivations which drive our 
interest in accurate low-energy tttt scattering experiments, 
eventually completed by independent experimental infor¬ 
mation of a different kind. In particular it is worth ex¬ 
plaining why a small difference in the values of scattering 
lengths like the one mentioned in the introduction may 
be relevant. The reader who is mostly interested in the 
analysis and less in the interpretation of the new data 
can skip this section and proceed directly to the following 
one. Nothing in our analysis depends on theoretical ideas 
summarized here. 

During the last two years our understanding of the 
pattern of chiral symmetry breaking in QCD has consid¬ 
erably evolved [11,12,13,14,15,16], in particular concern¬ 
ing the dynamical role of the number of light flavours. 
TTTT scattering provides important information about the 
SU{2) X SU(2) chiral structure of QCD vacuum in the 
limit of massless u and d quarks. The theoretical interpre¬ 
tation of this information, however, remains incomplete 
unless one learns how to disentangle the influence of the 
strange quark. If the strange quark were massless, we deal 
with SU{3) X SU{3) symmetry, which is the theoretician’s 
paradise [17]. The reason is that all remaining quarks (c, 
b, t) are heavy compared to the QCD scale, and conse¬ 
quently their influence on the chiral structure of the vac¬ 
uum remains tiny. For the same reason, we would reach 
another (a priori different) paradise if the strange quark 
mass were sent to infinity. In these two limiting cases, the 
question of the pattern of chiral symmetry breaking re¬ 
duces to the question of condensation of massless qq pairs 
in the vacuum; tttt scattering would by itself suffice to 
detect and fully control this one-fermion loop effect. Un¬ 
fortunately, the deplorable fact that we live in neither of 
these paradises [17] complicates the problem a bit. In the 
real world the strange quark is considerably heavier than 


u and d quarks and yet it is light compared to the QCD 
scale. Consequently, the vacuum gets polluted by massive 
virtual ss pairs. One may wonder how this fact could af¬ 
fect the SU{2) X SU{2) chiral structure of the vacuum, 
which is merely a matter of massless u and d quarks. In¬ 
deed, there would be no such effect in the large- A’c limit, 
in which the transition: 

ss < —> uu + dd ( 1 ) 

is forbidden and, consequently, strange and non-strange 
virtual pairs remain uncorrelated. Today we know that in 
the real world, the OZI-rule violating transition Eq. (1) is 
in fact rather important precisely in the vacuum channel 
[14] and that the resulting symmetry breaking correlation 
{ss{uu + dd)) affects the SU{2) x SU{2) chiral structure 
of the vacuum. We call this phenomenon which proceeds 
via at least two fermion loops the induced quark conden¬ 
sate [16]. It is proportional to iris, it enhances the effect 
of the genuine condensate of massless quarks characteris¬ 
tic of the ideal world with SU{3) x SU{3) symmetry, and 
it would persist even if the latter would be absent. The 
situation can be summarized by the formula: 

U{2) = S{3) TOs-^scalar j (2) 

where S{Nf) denotes the v.e.v. — {uu) of the lightest quark 
in the limit in which the first Nf quarks become massless, 
whereas the second term represents the induced conden¬ 
sate with ^scalar > 0 proportional to the amplitude of the 
transition Eq. (1) in the scalar channel. Since T’(3) is de¬ 
fined in the limit = md = nua = 0 and since there are 
no more massive quarks left which would be sufficiently 
light to pollute the vacuum, T’(3) = Eg Rg (in the current 
yPT notation) represents the true condensate of massless 
quarks of the SU{3) x SU{3) symmetric world. On the 
other hand, 27(2), defined at = md = 0 and the phys¬ 
ical value of ms , is detected in low energy tttt scattering. 
Our original quest for the importance of quark condensa¬ 
tion [18,19] did not sufficiently emphasize the subtlety of 
the theoretical and experimental distinction between 27(2) 
and 27(3). Phenomenological indications in favour of a 
substantial OZI-rule violating transition, Eq. (1), were not 
yet available and the large-Nc wisdom 27(2) ~ 27(3) dom¬ 
inated our thinking. Today, we have to answer the more 
precise question whether the suppression of the quark con¬ 
densate - expected in QCD due to the screening effect of 
multiflavour massless quark loops - is already visible for 
Nf = 3. 

The important consequence of the new E865 data [1] 
is that 27(2) is sufficiently large to keep the two-flavour 
COR ratio 27(2) = (m„ -f md)S{2)/F^M^ close to one 
[20] (see Sec. 6 for more details). However, we still do not 
know whether the observed size of 27(2) reflects an im¬ 
portant contribution of the genuine condensate 27(3) or 
whether it is (at least partially) due to the induced con¬ 
densate and to the strange quark mass, i.e. to the second 
term in Eq. (2). In the latter case, the genuine condensate 
27(3) could be rather small and the corresponding three- 
flavour COR ratio 27(3) = (to„ -I- md)27(3)/F^M.^ sig¬ 
nificantly below one. [Since 27(2) ~ 1, the suppression of 
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X (3) would imply an important contribution of rris to 
through terms mms,mm1 ..where m = {mu + md)/2] 
Moussallam [12] has proposed a rapidly convergent sum 
rule that allows an estimate of the transition Eq. (1) and 
of the size of the induced condensate. Detailed studies of 
this sum rule [12,13,15] are compatible with a suppres¬ 
sion of the condensate by a factor ~ 2 when going from 
2 to 3 flavours: ^(3) ~ Ar(2)/2 indicating that the two 
contributions to 27(2) on the right-hand side of Eq. (2) 
should be of comparable size. It is rather difficult to make 
a precise quantitative statement here: the evaluation of the 
Moussallam’s sum rule involves a few steps that include 
some unproved properties of the scalar form-factors and 
a certain model-dependence (of a similar type to the one 
in the phenomenological evaluation of the scalar radius of 
the pion [21]). Hence, it is crucial to find an independent 
way of disentangling the genuine and induced condensate 
contributions to 27(2). 

An analysis of tttt scattering exclusively based on SU (2) 
xSU{2) yPT can never separate the two components of 
27(2) in Eq. (2), for the simple reason that it does not 
involve any information on the actual size of m^. On the 
other hand, an SU{3) x SU{3) analysis of suitable tttt ob¬ 
servables supplemented by additional observables, such as 
Mk or Ffc, directly sensitive to mg, can considerably re¬ 
strict the possible values of three important parameters: i) 
the three-flavour GOR ratio X{3) = 2mE{3)/F^M^, ii) 
the quark mass ratio r = mg/m and Hi) the pseudoscalar 
decay constant Fq = EVlm„=md=ma=o- An example of a 
relation of this type is the strong correlation which exists 
between r and X{2) [11,15]. Further restrictions can be 
obtained from the emerging analysis oi K — tt scattering 
[ 22 ]. 

We are going to present this combined analysis in a 
separate paper [23]. It will be shown in particular that in 
order to be conclusive, one has to use extremely precise 
values of subthreshold tttt parameters a and l3 as the in¬ 
put of our analysis. Typically, the actual precision driven 
by the new E865 data analyzed in the following section 
will be just sufficient to reach conclusions about the sup¬ 
pression of the genuine condensate 2f(3) on the 1-a level. 
Notice that a more precise measurement of low-energy tttt 
scattering is conceivable in more dedicated experiments 
which are either ongoing [24] or planned [25]. 


3 The two S-wave scattering lengths 

3.1 Extended solutions of Roy equations 

The E865 data can be analyzed using a parametric rep¬ 
resentation of the solution of Roy equations for the tttt 
phase shifts. A set of dispersion relations derived by Roy 
[4] allows one to relate the phase shifts Sj), S) and S§ in 
the region 4M^ E s < sq (with yGo ~ 0.8 GeV) to data 
at intermediate energies (yso < \/s < 2 GeV) and to two 
subtraction constants. The latter can be identified with 
the two S-wave scattering lengths Oq and Qq. 

The Roy equations, analyzed thoroughly in Ref. [5], 
yield a boundary value problem. The solutions must in¬ 


terpolate between the phase shifts at the threshold, fixed 
by Oq and Oq, and the three phases at the matching point 
So: 


^0 —'^o(so)) ^1 —'^i(so)) ^2—< 59 ( 50 ), (3) 

determined from data above so. As stated in Ref. [5], the 
behaviour of the phases above so is less important than 
the boundary values, because they only affect the slope 
and the curvature of the solutions. With the experimental 
input encountered in practice, the system of Roy equations 
admits a unique solution provided that the matching point 
So is carefully chosen (0.78 GeV < ^/so < 0.86 GeV). 

Moreover, for given boundary conditions (og, 6 * 0 , 6 * 1 , 6 * 2 )) 
arbitrary values of ag generate a strong cusp in the P- 
wave solution at the matching point sg. If we require the 
phases to be smooth, Og is determined as a function of 
(og, 00 : ^ 2 ). Since 0 o,i ,2 can only vary in their experi¬ 

mental range, this requirement leads to a correlation be¬ 
tween Og and Og, defining the so-called Universal Band 
(UB) in the ( 09 , 09 ) plane. Different choices for 0 o.i ,2 rep¬ 
resent lines in the UB, Og = E(og). Inverting the rela¬ 
tion between Og and (og, 0 g, 0 i, 02 ), we can consider 02 
as a function of the other parameters. This means the 
solutions of the Roy equations do depend eventually on 
(ag,Oo, 0 o, 0 i) only. 

The data in the / = 0,1 channels lead to: 

0g = 82.3° ± 3.4°, 01 = 108.9° ± 2°. (4) 

The authors of Ref. [5] have provided explicit numerical 
solutions of the Roy equations for 0g = 82.0° and 0i = 
108.9°. We have included in the parametrization of Roy 
equations’ solutions^, an explicit dependence on 0g and 
01 , generating solutions (with the same driving terms and 
experimental input above the matching point sg) for nine 
different sets (0o,0i) and a few tens of (ag,ag) inside the 
UB. 

Following Ref. [5], we parametrize our solutions, for 
energies below 800 MeV, as: 


tan dj (s) = \ I — 


4M2 


-q 


\ s-si 


+ Cl 


(5) 


The dependence on Og and Og of the Schenk parameters 
X = A,B,C,D in Eq. (5) is well approximated by: 


Xi = Zi + Z2U -f Z3V + ZiU^ + z^v^ -\- ZqUV 

+zru^ + ZgU^V + Zguv'^ + ZioV^ , ( 6 ) 

where: 

qP qP 

u=-2-l, t,= -0-l, pg = 0.225, P 2 =-0.03706, 

Po P2 

_ (7) 

^ We have used the standard routines of MINUIT for all the 
minimization and fitting procedures of this paper. 
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while, for each coefficient Zi, the dependence on the phase 
shifts at the matching point is parametrized by: 


Zj = Gj + SOq bj -|- 59i c 


'3 > 


where 


500 = 6*0 - 82.3° , 501 = 0i- 108.9° 


( 8 ) 


(9) 


We first perform a fit using the solutions of the Roy 
equations of Ref. [5]. The is defined as: 

9 /(^2)ACGL(^exp)_(^2)expy 


exp 


Xgiobai(ao,ao) = 

j=l \ '^3 

“ ^[ 50 -< 5 l]ACGL(,exp)_j^O_^l]exp) 


E 


exp 


( 11 ) 


The parameters Sg, s} and Sg are fixed by the boundary 
conditions: 

'^o(so) = ^0 7 '^l('5o) = Si , 1^0(^o) = 6*2 , (10) 

where 6 l 2 (ag, Ug, 6 lg, ^i) is parametrized following Eqs. ( 6 ) 
and (8). The coefficients aj, bj, cj are collected in App. B. 
This parametrization describes our solutions to better than 
0.3° for the I = 0,2, and 0.5° for the / = 1 partial waves 
in the Universal band. 

By setting 6 *g = 82.0° and 0i = 108.9°, we can com¬ 
pare with Ref. [5]. We obtain slightly different Schenk pa¬ 
rameters for the so-called reference point Ug = 0.225 and 
Og = —0.0371, but the phase shifts are identical up to 
a few tenths of a degree. We obtain the same Universal 
Band, and only its lower half meets the consistency con¬ 
dition (Roy equations fulfilled in their range of validity 
2M,r < y/s < ^/si = 1.15 GeV). In the range of interest 
for Og, the gap between the parametrization of Ref. [5] and 
our Roy solutions is at most 0.3° in the / = 0,2 channels 
and 0.7° in the 1 = 1 channel (for -y/s ~ 0.7 GeV and 
much smaller near threshold in all the channels). 


3.2 Model-independent determination of Ug and Ug 

The E865 data on Sq — SI were analyzed in Ref. [1] in order 
to extract Ug. While two different (although compatible) 
results were quoted in this reference for Og, (og = 0.228 ± 
0.012±0.004l°;°'j’® and a[j = 0.216±0.013±0.004±0.005), 
no results were given for Og, which is harder to pin down 
from Ke 4 data alone (see Sec. 1). An unconstrained fit of 
E865 data, using the Roy equations of Ref. [5], leads to 
a rather inaccurate result for the 1 = 2 scattering length: 
ag = 0.237 ± 0.033 and a§ = -0.0305 ± 0.0226 (y^ = 
5.44/5 d.o.f and the correlation coefficient is 0.96). 

In order to extract both scattering lengths additional 
information has to be provided. The first constraint arises 
from the necessary consistency of the Roy solutions with 
the 1 = 2 data above the matching point. This forces the 
S-wave scattering lengths to lie within the so-called Uni¬ 
versal Band. Unfortunately, this model-independent con¬ 
straint is rather weak. 

We make use of additional information by Htting a 
broader set of data below 800 MeV, namely Rosselet and 
E865 sets for 7 = 0,1 [2,1] and Hoogland (sol. A) and 
Losty sets for 7 = 2 [6,7]. Notice that a similar fit has 
been considered in Refs. [5,26] but without the E865 data 
(cf. Figs. 11 and 12 in Ref. [5]). 


where [(5/]^®®^(aQ, Og, s) is the parametrization of Roy so¬ 
lutions proposed in Ref. [5]. i and j are the indices of the 
experimental points for 7 = (0,1) and 7 = 2 respectively. 
This fitting procedure, referred to as “global”, yields: 

Global: a]] = 0.228 ± 0.012, al = -0.0382 ± 0.0038. 

( 12 ) 

with Xmin = 16.45/18 d.o.f and a correlation coefficient 
of 0.788. Including data on the P-wave below 800 MeV 
reported by Protopopescu et al. [27] does not change the 
result of the fit, yielding^ y^j^^ = 23.1/28 d.o.f. and Og = 
0.228 ± 0.012, al = -0.0392 ± 0.0038. 

A second fitting procedure can be followed, in which 
we use our solutions of the Roy equations to include the 
dependence on the phase shifts at the matching point 
(0g, 0i). The y^ is then dehned as: 


9 

Xextended(® 0 ’ ^ 0 ? ^l) ^ ( 




exp 

i _ 


2 


+ 



-82.3°y 

' 3 ( 4 ° ) 





- 108.9° \ 
J 


2 


(13) 


where [(5/]°’°*(ag, Og, 0o, 0i,s) is our extended parametriza¬ 
tion of Roy solutions. This Ht, called “extended”, leads to 
the same S-wave scattering lengths as the “global” fit: 

Extended: ag = 0.228 ± 0.013, al = -0.0380 ± 0.0044, 
6 ig = 82.1° ±3.3°, 6»i = 108.9° ± 2.0° 

(14) 

with Xmin — 16.48/18 d.o.f and the correlation matrix: 



Og al 00 

01 

al 

1.000 0.799 -0.319 

-0.004 

al 

- 1.000 -0.271 

0.029 

00 

- 1.000 

0.000 

01 

_ _ _ 

1.000 


The results of these analyses are shown in Fig. 1, where 
we have indicated the 1- and 2-cr contours for both de¬ 
terminations. These contours are defined respectively as^ 

^ Specifically we have used the energy-independent solution, 
which has larger error bars. The P-wave production data at 
low energy have been often criticized, most recently in Ref. [5]. 
Therefore we prefer not to include them in the analysis, but 
instead use our results to predict the P-wave at low energy. 

® We recall that in the case of the simultaneous determina¬ 
tion of two variables, this definition of the contours correspond 
to 39% and 86% confidence level. 
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= Xmin + 1 and = Xmin + 4. We see that our fitting 
result lies slightly below the center of the UB, where the 
consistency condition for Roy equations is met. 

A comment is in order here concerning our use of the 
ACM (A) 1 = 2 phase shifts by Hoogland et al. These were 
extracted following the method A, which is a conventional 
Chew-Low extrapolation to the pion pole of the measured 
f-channel (m = 0) helicity moments [6] (the beam momen¬ 
tum was 14 GeV/c). A similar method was used by Losty 
et al. [7]. A second method (B) is presented in Ref. [6], 
which is based on an overall ht of an (absorption) model 
for the amplitude to all non-negligible s-channel helicity 
moments (m = 0 and m = 1). The method B involves 
extra assumptions and parameters, some of which exhibit 
unexpectedly rapid energy variations. No x^ is given in 
Ref. [6] (in a preliminary analysis [28], based on a third of 
the data, a poor x^ was reported for method B). 

We have tried to use in our fit solution B of Ref. [6] 
instead of solution A. Due to the small error bars of the 
former, we did not succeed in obtaining a consistent de¬ 
scription of both ACM(B) and E865-Are4 data within the 
ACGL solutions of Roy equations [5]. The minimum has 
X^ = 68/18 d.o.f. and is situated far outside the Univer¬ 
sal Band. Such a Ht has little meaning, since the AGGL 
solutions are valid exclusively inside the Universal Band. 
Solution ACM(A) is free from such difficulties. 

It has been suggested in Ref. [5] that the difference 
between the phase shifts AGM(A) and (B) indicates size¬ 
able systematic errors, and that the errors associated with 
ACM(A) solution should consequently be enlarged. It is 
not obvious to us that method B yields a correct estimate 
of the systematic errors inherent to method A - especially 
since the two methods do not use the same sample of data. 
We find it nevertheless useful to show in App. A and in 
Fig. 6 how our results would be modified if the errors in 
ACM(A) phase shifts were increased according to the pre¬ 
scription advocated in Ref. [5]. Let us mention here that 
these modifications barely affect the main conclusion. 

3.3 Comparison with the xPT prediction for 2ag — bog 

In the theoretical prediction of Og and Og based on stan¬ 
dard xPT including 0(p®) accuracy [9,10,20,29], one may 
distinguish two steps. The first step concerns the relation 
between the combination 2ag — bog and the scalar radius of 
the pion (r^)s [10,29]. This step is practically independent 
of the badly known 0 {p‘^) constant £3 but it requires an in¬ 
dependent phenomenological determination of {r'^)s and it 
is rather sensitive to the two-loop corrections (a more de¬ 
tailed discussion of this theoretical prediction can be found 
in Sec. 7). If one takes the value {r^)s = 0.61 ± 0.04 fm^, 
the prediction amounts to a narrow strip in the Og — Og 
plane, given in Ref. [10] and reproduced in Fig. 1: 

al = G(a[;) ± .0008, (16) 

where the function G(ag) may be parametrized as 

G{al) = -.0444 -k .236(a[; - .22) - .61(a[; - .22)^ 

-9.9(a°-.22)3, (17) 


and the error bar is estimated within the theoretical frame¬ 
work defined in Refs. [10,29] (see Sec. 7 of the present 
paper). The second step of the prediction procedure then 
consists in locating the actual position inside the narrow 
strip Fq. (16) and it involves, among other things, an es¬ 
timate of the constant £ 3 . 

The analysis performed in the previous subsection ma¬ 
kes use only of the F 86 b and Rosselet data on (5° — <5^, the 
Hoogland and Losty data on $1 together with the solution 
of Roy equations, and does not use xPT or a particular 
value of {r‘^)s- It provides thus a sensitive experimental 
test of the theoretical prediction represented by the CGL 
correlation Fqs. (16) and (17). It is seen from Fig. 1 that 
the l-cr ellipses resulting from both fits “global” and “ex¬ 
tended” are situated outside the CGL narrow strip despite 
the fact that they are entirely contained within the Uni¬ 
versal Band required by the consistency of Roy equations 
solution. On the other hand, the 2-cr contours intersect the 
narrow strip. We thus conclude that there is a marginal 
(l-cr deviation) discrepancy between the theoretical pre¬ 
diction Eq. (17) and the result of the “global” and “ex¬ 
tended” fits. This discrepancy will be further commented 
on in Sec. 7. 

In order to make the comparison more quantitative, 
we can perform a fit to the E 86 b data alone, imposing 
by hand the correlation described by the narrow strip. A 
similar fit has been performed in Ref. [10], leading to the 
central value ajj = 0.218 (no uncertainty was indicated in 
that reference). Our fitting procedure is defined by: 



where k runs only through the E 86 b points. We have ob¬ 
tained: 

Scalar: a[] = 0.218 ± 0.013, al = -0.0449 ± 0.0033. 

(19) 

with Xmin = b.89/b d.o.f., compatible with the results of 
E 86 b [1] for Og, and a correlation coefficient of 0.972. The 
corresponding l-cr contour is indicated in Fig. 1. We refer 
to this fit and to the corresponding l-cr ellipse as “scalar” 
to be compared with the “global” and “extended” fits. 
The meaning of x^ and of the standard deviation in the 
“scalar” fit should be taken with caution: The error bar 
0.0008 arises from uncertainties in the experimental in¬ 
put, while the theoretical errors inherent in the estimate 
of 0(p®) corrections are more difficult to quantify. On 
the other hand, the fits “global” and “extended” are fully 
based on experimental data and corresponding errors. 

Finally, we would like to briefly comment on the Olsson 
sum rule for 2aQ — Sog, as discussed in Ref. [5]. This sum 
rule converges slowly and demands good control over the 
asymptotic contribution, which is hard to obtain outside 
specific models. According to the model used for this pur¬ 
pose in Ref. [5], the asymptotic contribution to 2aQ — Sog 
is Gas = 0.102 ± 0.017. Even with such small error bar, 
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Fig. 1. Fit results using either the Roy solutions of Ref. [5] (red, dotted ellipse - “global”) or our parameterization of the Roy 
solutions (black, solid ellipse - “extended”). In each case, the thicker lines indicate the 1-a ellipse, and the thinner ones the 2-cr 
contour. The Universal Band (delimited by the two straight lines) is drawn according to Ref. [5], and the narrow strip (shaded 
region, cyan) related to the scalar radius of the pion is taken from Ref. [10]. We have indicated the result of the fit “scalar” 
(blue, hatched ellipse) and the SxPT prediction of Ref. [29] (red, filled small ellipse). 


the final result shown in Eq. (11.2) of Ref. [5] is consis¬ 
tent with our global fit, which leads to ( 2 oo — 5oo)giobai = 
0.647 ±0.015. If the actual error bar in Oas is bigger, the 
impact of the Olsson sum rule on our fit becomes com¬ 
pletely negligible. 


4 Subthreshold parameters 


In the low-energy domain the tttt amplitude is strongly 
constrained by chiral symmetry, crossing and unitarity. As 
was first shown in Ref. [18], the amplitude depends on only 
six parameters up to and including terms of order [jp/A h)^ 
in the low-energy expansion. In Ref. [ 8 ], the amplitude was 
written as: 


AKMSF(s|t, u) = A™*(s|t, u) ± u) 


( 20 ) 


AP°\s\t,u) = 


A 

F3 


4M.f 


s — 


FI 3 


±4(s - + 2Mlf + {u- 2Mlf] 


pA' 

-*■ TT 

^3 

Fr 


±4 

-*■ TT 

A 

FB 


± + 2Mlf ± (u - 2Mlf] 


(21) 


where is a known function of the Mandelstam vari¬ 
ables s,t,u that collects the unitarity cuts of the ampli¬ 
tude and explicitly depends on a, (3, Ai, A 2 . 

This amplitude was constructed as the general solu¬ 
tion of unitarity, analyticity and crossing symmetry up to 
and including OijB). The six parameters a,/3, Ai,.. ., A 4 
correspond to an expansion of the amplitude in the cen¬ 
tral region of the Mandelstam triangle, and are therefore 
called “subthreshold parameters”. A complete calculation 
in the framework of Standard yPT [9] confirmed this re¬ 
sult, allowing in addition to relate the six parameters to 
the quark masses and EEC’s of the standard chiral La- 
grangian. This last step is crucial to translate the exper¬ 
imental information into knowledge of the EEC’s, which 
parametrize the chiral structure of the vacuum of QCD. 
The six parameters introduced in Ref. [9], 61 ,... , 65 , are 
dimensionless combinations of EEC’s in one-to-one (lin¬ 
ear) correspondence with a,/3, Ai,..., A 4 of Ref. [8], or 
with Cl,..., Ce, subsequently introduced in Ref. [20]. Dif¬ 
ferent choices for the set of six subthreshold parameters 
correspond to different parametrizations of solutions of 
unitarity, analyticity and crossing symmetry constraints, 
which are equivalent up to 0 {p^) and only differ at 0 {p^). 

On the other hand, the Roy equations allow one to 
determine the low-energy amplitude in terms of only two 
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subtraction constants, identified with the two scalar scat¬ 
tering lengths. It is therefore possible to match the two 
amplitudes in their common domain of validity, in order 
to determine, through the experimental determination of 
the scattering lengths, the six subthreshold parameters. 
Such a program was already advocated in Ref. [30], lead¬ 
ing to rapidly convergent sum rules for the parameters 
Ai,..., A 4 . A similar matching procedure, using new solu¬ 
tions of the Roy equations, has been carried out in Ref. [29] 
(subtracting the dispersion integrals at s = 0). Let us 
briefly outline the various steps, using the notation and 
results of this last reference. 

Starting with particular values of Og and Og, we can 
use the solutions of the Roy equations to compute the low- 
energy moments J^. In conjunction with the background 
moments and H, estimated in Ref. [5], we compute the 
phenomenological moments and their linear combina¬ 
tions denoted .6 and defined in Eq. (3.5) of Ref. [29]. 
Matching the phenomenological and the chiral representa¬ 
tions of the amplitude connects the phenomenological pa¬ 
rameters Pi=i ,..6 and the chiral ones Ci=i ,..6 [see Eq. (4.2) 
of the same reference]. App. A and B of Ref. [29] can then 
be used to translate the chiral parameters Ci=i ,..6 into the 
parameters 6 i=i ...6 defined in Ref. [9], and finally into a, 
P and Ai=i,.. 4 . 

We can repeat this procedure for each set of (ug, Ug), as 
determined from the “global” and “scalar” fits, or (ao,aQ, 
0o,9i), as determined from the “extended” fit, described 
in the previous section. In order to take full account of 
the theoretical and experimental correlations among the 
six parameters, we proceed in the following way: we gen¬ 
erate random sets of (og, Og) or (ug, Og, 60 , 9i), distributed 
according to the 2- or 4-dimensional gaussian obtained 
from the covariance matrix of the fit. We then fit the re¬ 
sulting distributions for the subthreshold parameters by 
gaussians, leading"^ to Table 1. pa/s denotes the correla¬ 
tion coefficient between a and /3. 

The slightly larger error bars of the “extended” fit, 
compared to the ones of the “global” fit, reflect the in¬ 
fluence of the uncertainties in 0 g and 9i, which in the 
“global” fit are not explicitly taken into account. The dif¬ 
ferences in the central values between these two fits (al¬ 
though compatible within the errors) may be ascribed to 
the fact that the “extended” parametrization is not as ac¬ 
curate as the ACGL one, due to the fact that the former 
has to account for the dependence on two more variables. 
The column referring to the “scalar” fit should be un¬ 
derstood as originating from a mixture of E865 data and 
yPT-based theoretical predictions that rely on assump¬ 
tions about the size of 0(p®) counterterms (see Sec. 7). 
For this reason the associated errors should not be inter- 


^ The phenomenological moments are integrals of the 
I = 0,1,2 phase shifts from threshold to y/s 2 = 2 GeV. The 
solutions of the Roy equations are used for s < so, and experi¬ 
mental input is used above the KK threshold. An interpolation 
is necessary in the intermediate region [so,4M^]. We have ob¬ 
served a weak sensitivity of Ai and A2 on the interpolation 
prescription. On the other hand, the values of a and /3 are 
independent of this procedure. 


Table 1. Subthreshold parameters for the three different fits 
considered in this paper. See the text for a discussion of the 
error bars. 



Global 

Extended 

Scalar 

61 

-2.08 ± 6.12 

-1.51 ± 7.01 

-13.06 ± 5.31 

&2 

9.35 ± 1.43 

8.93 ± 1.62 

11.51 ± 0.62 

63 

-0.38 ± 0.03 

-0.36 ± 0.07 

-0.33 ± 0.01 

64 

0.716 ± 0.008 

0.710 ± 0.010 

0.727 ± 0.005 

65 

3.21 ± 0.25 

3.21 ± 0.44 

3.53 ± 0.15 

be 

2.23 ± 0.07 

2.20 ± 0.08 

2.34 ± 0.03 

a 

1.381 ± 0.242 

1.384 ± 0.267 

1.034 ± 0.248 

p 

1.081 ± 0.023 

1.077 ± 0.025 

1.116 ± 0.010 

Paf 3 

-0.14 

-0.23 

0.53 

Ai ■ 10® 

-4.40 ± 0.28 

-4.18 ± 0.63 

-3.97 ± 0.12 

A2 • 10® 

9.04 ± 0.10 

8.96 ± 0.12 

9.17 ± 0.06 

As • 10^ 

2.21 ± 0.10 

2.22 ± 0.16 

2.32 ± 0.06 

A4 • 10^ 

-1.40 ± 0.04 

-1.38 ± 0.04 

-1.46 ± 0.02 


preted in the strict statistical sense. The corresponding I- 
and 2-cr ellipses in the (a — /?) plane are drawn in Fig. 2. 

It is worth stressing that a rather small difference be¬ 
tween Ug and flg resulting from the scalar fit, on the one 
hand, and from the global and extended fits on the other 
hand, results in a more pronounced difference in the cor¬ 
responding values of the subthreshold parameters a,p. 
Whereas the scalar fit (and the CGL prediction) is char¬ 
acterized by values of a close to (or smaller than) I and 
P well above l.IO, the global and extended fits lead to 
central values of a ^ 1.4 and relatively smaller values of 
p. It will be shown elsewhere [23] that this qualitative 
difference finds its interpretation within the three-flavour 
analysis of tttt scattering together with other observables. 



Fig. 2. Gorrelation between a and P from the fit “scalar” 
(blue, dashed ellipse) and our fits (black, solid for the “ex¬ 
tended” and red, dotted for the “global”). The thicker lines 
correspond to the l-cr ellipses, whereas the thinner ones indi¬ 
cate the 2-0 ellipses (not shown for the “scalar” fit). 
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5 Nf = 2 Mass and Decay constant Identities 

In order to investigate the consequences of the results 
obtained so far for the parameters of the effective La- 
grangian, we start with the Ward identity satished by the 
two-point function of the divergence of the axial current 
uifilsd and of its conjugate at zero momentum transfer. 
We isolate all LO (linear) and all NLO (quadratic) contri¬ 
butions in the quark mass m = (to„ -I- md)/‘2,, and collect 
all 0{m^) and higher order terms into the NNLO remain¬ 
der d: 


F^Ml = + ^(44 - 4 ) + F^M^S. ( 22 ) 

Here 

M^ = 2mB, B = IJ{2)IF^ (23) 

are defined in the SU{2) x SU{2) chiral limit, keeping the 
strange quark mass at its physical value: 

F = lim .^TTI physical 5 (^4) 

^(2) = lim (w-'w) 1777,^—physical ■ (25) 

rriu 


A similar identity holds for the two-point function of axial 
minus vector currents, giving 


F^ = F^ + ^h + F^e. (26) 

On the right hand side, all LO and NLO contributions 
are again explicitly shown, and all higher orders - 0(rn?) 
and higher - are included in the NNLO remainder e. £3 
and £4 are the standard SU{2) x SU(2) LEC’s, which are 
scale-independent and exhibit logarithmic singularity in 
the chiral limit. They can be defined non-perturbatively, 
via the low-energy behaviour of the two-point functions 
that enter the Ward identities considered above. It may 
be useful to consider the NNLO remainders d and e in 
Eqs. (22) and (26) as known, and to treat the mass and 
decay constant identities exactly, not performing any ex¬ 
pansion. This avoids the use of perturbation theory when 
eliminating the order parameters and F^ in favour of 
observable quantities. We will see later that, despite the 
fact that q: — 1 need not be particularly small, this non- 
perturbative precaution is not absolutely necessary in the 
Nf = 2 case. It will however, be fully justified in the case 
of three light flavours. In any case, the above method - 
which does not coincide either with the Standard or with 
the Generalized yPT - is completely meaningful no mat¬ 
ter how large £3 or how small the condensate might be. 

The fundamental order parameters, in appropriate units, 
the condensate and the decay constant defined as: 


Xi 2 ) 


2mS{2) 

~FfMf 


Y{2) 


2mB 

~mF 



X{ 2 ) 

Y{2) ’ 
(27) 


are related to the observables Mt^,Ft^, to the LEC’s £3 
and £4 and to the NNLO remainders S, e by the following 


identities: 

w.oj = _ 2(1-J) _ 

^ ’ 1 -e+[(l-e) 2 -24^(1 - ’ 

X{2) = 1-6- (44 - 4)?T(2)V2 , (29) 

Z{2) = l-s- 2£4(Y{2) . (30) 


In the last three equations, we have denoted 

1 

^ 167r2 F2 ■ 


(31) 


The TTTT subthreshold parameters a and (3 can be expressed 
similarly. Reading the LO and NLO perturbative contri¬ 
butions to F^M^a and to F^/S from Ref. [31], one obtains 
the identities: 


a = 1 - (1 -P 34 - 44)CE(2)2/2 -k , (32) 

P=l + 2(4 - 1)C4"(2) + £/3. (33) 

It may be expected - at least for the moment - that the 
NNLO direct remainders 6 a and s /3 are not more impor¬ 
tant than the uncertainties in the determination of the 
parameters a. and /3 from experimental data. 


6 Determination of LEC's and order 
parameters 

If we expand the previous expressions of a and (3 in powers 
of 4 we obtain the following (linearized) expressions in 
term of the LEC’s 4 and 4: 

a-/3 = 34l-4)/2, (34) 

£3-1 = 244 - 1). (35) 

This is an excellent approximation, unless £3 or £4 become 
“too” large. However, even if one of them were large, the 
non-linear equations (32) and (33) of the previous section 
would still be exact identities; moreover, the definition of 
4 and £4 in terms of two-point functions is independent 
of their magnitudes. We can use Eqs. (34) and (35) to 
translate our determination of (a, (3) into a l-cr contour 
plot in the £3 — £4 plane. In Fig. 3, we show three ellipses 
corresponding to those in the a — (3 plane displayed in 
Fig. 2 above. 

If we use the formulae Eqs. (28), (32) and (33), but now 
without linearizing, the previous ellipses are deformed, as 
shown in Fig. 4 (solid lines). The corresponding contours 
in the X{2) — Z{2) plane are plotted in Fig. 5. Up to 
now, we have neglected the indirect remainders 6 and e 
as well as direct remainders 6 a and e^. In the case Nf = 
2, we expect these NNLO quantities to be less than 1%, 
since Eqs. (22) and (26) are obtained by an expansion 
in powers of the nonstrange quark mass. This leads to a 
(small) additional broadening of the l-cr regions, as seen 
in the plots with thinner lines in Figs. 4 and 5 (do- and ep 
are negligible compared to the present uncertainty in the 
parameters a and (3). 
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It is clear from Fig. 4 that we obtain rather large and 
negative values of £ 3 , compared to the standard expecta¬ 
tion of 2.9 ± 2.4. This can be interpreted as a consequence 
of an important OZI-rule violating transition Eq. (1) lead¬ 
ing to a larger value of the Nf = 3, large-lVc suppressed 
constant Lq{Mp), than usually assumed. Sum rule esti¬ 
mates of the latter [12,15] support this interpretation. 
Notice that the OZI rule is an important ingredient of 
the standard estimates of I 3 [3]. We will return to this 
question elsewhere [23]. 

In Fig. 5, we see that the two-flavour GOR ratio ^"(2) 
is constrained (at one sigma): X{2) = 0.81 ± 0.09. In 
Ref. [32], this ratio has been theoretically analyzed, in¬ 
cluding the 0(p®) double chiral logarithms of Generalized 
yPT, determining X{2) as a function of a -I- 2/3 and other 
leg’s, a combination of the results of the latter analysis 
with the present (correlated) values for a and /3 leads to a 
range of values for X{2) completely consistent with ours. 
Further examination of Fig. 5 shows that is also 

limited to a rather narrow band, (F/F^.)^ = 0.90 ± 0.03. 
Let us mention that the difference seen in the (a, /3) plane 
between the various fits reappears clearly here. The scalar 
fit (and the GGL prediction) favours larger values of X(2) 
and lower values of Z{2) than the global/extended fit. 



3 

Fig. 4. l-a ellipse in the £3 — £4 plane, from the fit “scalar” 
(blue, dashed ellipse) and our fits (black, solid for the “ex¬ 
tended” and red, dotted for the “global”), using the non- 
linearized formulae. The thinner lines indicate the domains 
allowed if e, 5 < 1%. 



-3 

Fig. 3. l-a ellipse in the £3 — £4 plane, from the fit “scalar” 
(blue, dashed ellipse) and our fits (black, solid for the “ex¬ 
tended” and red, dotted for the “global”), using the linearized 
formulae. 


7 Comments on the correlation between 
scattering lengths and the scalar radius of 
the pion 

A relation between 2aQ — fiog and the scalar radius of the 
pion based on two-loop yPT has been derived in Refs. [10, 
29]. For the current value of the scalar radius (r^)s = 
0.61 ± 0.04 fm^, this prediction results in the narrow strip 



Xf2j 

Fig. 5. l-a ellipse in the X{2) — Z{2) plane, from the fit 
“scalar” (blue, dashed ellipse) and our fits (black, solid for 
the “extended” and red, dotted for the “global”). The thinner 
lines indicate the domains allowed if e, 5 < 1%. 


in the Og — Og plane shown in Ref. [10] and reproduced 
here in Fig. 1. The accuracy of this prediction is not only 
conditional on the experimental error, but also on theoret¬ 
ical assumptions and “rules”, which are a priori reasonable 
and natural, but not necessarily true. 

First of all, these assumptions concern the estimates 
of the 0(p®) corrections [20,29]. Although it is gener¬ 
ally admitted that for certain observables (e.g. scatter¬ 
ing lengths), these corrections could be enhanced, it is 
expected that for quantities receiving a weak contribu¬ 
tion from O(p^) loops, the whole chiral expansion will 
rapidly converge. An example of this rule of thumb which 
is relevant for the present discussion is the quantity Ci — 
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M^(r^)s/3, which does not contain 0 {p‘^) logarithms at 
all. (Cl stands for a combination of subthreshold parame¬ 
ters defined in Ref. [29].) We shall argue that this fact does 
not prevent the corresponding 0(p®) corrections from be¬ 
ing relatively important. (An alternative “rule” is conceiv¬ 
able [23]: expand QCD correlation functions at kinemat- 
ical points sufficiently distant from all Goldstone boson 
singularities. Such rule is a priori not less or more natu¬ 
ral.) Next, it is usually assumed that 0{p^) counterterms 
at a suitable scale can be estimated via the narrow reso¬ 
nance saturation [9,29,33]. In fact, already at 0{p*) this 
assumption fails in channels where 1 /Nc - corrections are 
large and/or the OZI-rule is strongly violated. This is what 
likely happens in the scalar channel, which is particularly 
relevant for the present discussion. Furthermore, the exist¬ 
ing resonance estimates of 0(p®) counterterms have been 
so far based on a “resonance effective Lagrangian £res” in¬ 
volving (and missing) the same resonances with the same 
“minimal resonance couplings” as in Ref. [34], originally 
used to estimate the O(p^) LEC’s. It has often been ar¬ 
gued [35,36] that additional non-minimal couplings are 
necessary to avoid conflicts with the QCD short-distance 
behaviour of two- and three-point functions, especially if 
the latter involve (pseudo)scalar currents. The estimates 
of the corresponding 0(p®) counterterms ri, r 2 , ra, r 4 , rs 
can be affected by these new resonance couplings. We shall 
return to the resonance estimates of the r„’s shortly. 



Fig. 6. Correlation Eq. (36) between aj] and Uq, obtained 
from the relation between 2ao — 5ao and the scalar radius of 
the pion. The shaded oblique band (yellow) corresponds to 
the 0(p'*) level [Eq. (36) with Sa =0], using the dispersive 
estimate of Ref. [21]: (r^)s = 0.61 ± 0.04 fm^. The shaded 
curved strip (cyan) represents the effect of 0(p®) corrections 
Sa according to the prediction of Ref. [10] using the same input 
value of {r^)s- The contours correspond to fits using either the 
Roy solutions of Ref. [5] (red, dotted ellipse - “global”) or 
our parameterization of the Roy solutions (black, solid ellipse 
- “extended”), with the errors of ACM(A) data enlarged as 
detailed in Appendix A. In each case, the thicker lines indicate 
the 1-(T ellipse, and the thinner ones the 2 -(t contour. The other 
elements are identical to Eig. 1. 


Table 2. Values of Sa required for satisfying Eq. (36) and the 
estimate {r^)a = 0.61 ± 0.04 fm^. The correlation between uq 
and flo has been taken into account to compute the error on 
2ag -5ag. 



Scalar 

Global 

Oo 

0.218 

± 

0.013 

0.2279 

± 

0.012 

a-o 

-0.0449 

± 

0.0033 

-0.0382 

± 

0.0038 

2ag - 5ag 

0.660 

± 

0.011 

0.647 

± 

0.015 

Sa 

0.033 

± 

0.012 

0.020 

± 

0.015 


Finally a remark should be made about the dispersive 
estimate of {r‘^)s in Ref. [21], and the uncertainty related 
to it. The pion scalar form-factor and radius are not exper¬ 
imentally measurable quantities - information about them 
can only come from indirect theoretical constructions. In 
contrast to the case of the (observable) vector form-factor, 
QCD does not restrict very much the high momentum be¬ 
haviour of the scalar form-factor. It even does not guar¬ 
antee that the latter satisfies an unsubtracted dispersion 
relation. Consequently, the dispersive evaluation of (r^)s 
suffers from a certain model dependence concerning the 
higher momentum contributions; this is usually not dis¬ 
cussed in the literature. The quoted uncertainty should 
not be interpreted outside the framework of the model 
used in evaluating the scalar radius. 

Most of these critical remarks are obviously not new. 
It might however be useful to keep them in mind when dis¬ 
cussing the origin of the discrepancy between our model- 
free determination of scattering lengths from the data 
and the CCL narrow strip prediction [10]. The origin of 
the narrow strip Eq. (17) is more easily understood from 
Eq. (3) of Ref. [10]: 

20 ° - Soq = 

= 0.57158 -b 0.05541 (+ 5a ,(37) 

VO.GlfmV 

where 5a = 0{m^). It is worth stressing that the 0{p^) 
contribution 5a is essential for the numerical coherence of 
Eq. (36). If the 0(p®) contribution 5a is dropped, Eq. (36) 
reduces to the 0{p‘^) low-energy theorem [31] relating 2aQ— 
Soq and {r^)s- This model-independent 0{p'^) relation is 
represented in Fig. 6 as a straight oblique band corre¬ 
sponding to the estimate of the scalar radius {r‘^)s = 
0.61 ±0.04 W. 

In order to reproduce the prediction of Ref. [29]: 2aQ — 
5aQ = 0.663 ± 0.006, either the scalar radius should be as 
large as 1.01 fm^, or the 0(p®) correction 5a should move 
the 0 {p‘^) straight oblique band up to the curved strip 
along the bottom of the Universal Band, reproduced in 
Fig. 6 from Ref. [10]. If we believe the estimate {r'^)s — 
0.61±0.04fm^, andifEq. (36) should hold inside the scalar 
or the global ellipse, 5a should take the values indicated 
in Table 2. The 0{p^) correction amounts thus to 5% in 
the scalar case, and 3% in the global one. 

This size of 0(p®) corrections is consistent with gen¬ 
eral expectations. On one hand, NNLO contributions to 
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“smooth” observables - typically, QCD correlation func¬ 
tions far from Goldstone boson singularities - are expected 
at 1% level (see the NNLO remainders 6 and e: in Sec. 5). 
On the other hand, we expect an enhancement of higher- 
order corrections to infrared-singular “threshold quanti¬ 
ties”, such as scattering lengths [20,29,30]. Let us em¬ 
phasize that, in order to reproduce the narrow strip of 
Ref. [10], the two-loop correction 6 a has to be known 
very precisely: this is no longer a matter of a model- 
independent low-energy theorem [31] nor of an accurate 
knowledge of {r'^)s- In particular, the CGL prediction [10] 
would be brought into agreement with the present deter¬ 
mination of Oq and Oq by the “global” ht described in 
Sec. 3, if the 0(p®) correction Sa were reduced by factor 
2, even if the scalar radius were not modified. 

To understand the structure of the 0(p®) correction 
6 a, we follow Ref. [29] and write: 

2«r5«0 = ^C'i+M>i+O(MS). (38) 


Qfi is a combination of phenomenological moments defined 
in Eq. (6.4) of Ref. [29]. Ci collects polynomial coefficients 
in the chiral representation of the scattering amplitude. Its 
expansion in powers of m reads [29]: 


C, = l + ^{r^), + ^+eAi+0{e), (39) 

where 


Zii = 


-71L 



80; 


- + 44 + 


5393 


315 


252 3024 

contains the O(p^) constants 


3118- 79- 144; 

-I- 4 + 44 - 44 - 2fs2 ) 


(40) 


e„=£„-L, L = log ^, (41) 

and the scale-dependent 0(p®) counterterms 

4 = (47r)'‘r„(^). (42) 

Comparing with Eq. (36), the 0(p®) correction 6 a is 
seen to consist of two parts: 


6 a — 6 m + 1^1 • (43) 

6 m is essentially the phenomenological moment M^ai from 
which the leading infrared singularity has been subtracted: 

353 

6 m = M^ai - 47’‘—. (44) 

With the help of Eq. (6.8) of Ref. [29], one easily checks 
that in the chiral limit m ^ 0 all O(p^) contributions on 
the right-hand side of Eq. (44) exactly cancel. In practice, 


however, M^ai is obtained by evaluating infrared singular 
sum rules using physical scattering lengths and not their 
values in the chiral limit. Consequently, the quantity 6 m 
is expected to be sensitive to the infrared enhancement 
of higher order corrections to the scattering lengths. One 
obtains: 

M^il 

scalar — 0.0604±0.0054, (5 m [scalar = 0.034±0.0054, 

(45) 

M^ai [global = 0.0636 ± 0.005, (5 m [global = 0.037 ± 0.005, 

(46) 

for Oq and Oq inside the “scalar” and “global” ellipses re¬ 
spectively. The second part of the 0{p^) correction 6 a is 
due to Z\i: 

(5i = . (47) 

As shown in Eq. ^40), the latter depends on the O(p^) con¬ 
stants 1 ^ 1,^25 4)4 and on the counterterms 4 ) 4 ) 4)42 
describing 0 {p^) symmetry breaking effects in the scalar 
sector. The CGL prediction [10,29] of the strong correla¬ 
tion Eq. (17) between scattering lengths originates from 
a particular matching procedure worked out in Ref. [20], 
and from the resonance estimate of 0(p®) counterterms [9, 
29,33], whose practical outcome is the fact that (5i is neg¬ 
ligibly small and the whole two-loop correction 6 a reduces 
to the contribution 6 m as shown in Eq. (45). This conclu¬ 
sion is independent of the actual error in the dispersive 
estimate of the scalar radius {r‘^)s = 0.61 ± 0.04 fm^. 

The smallness of (5i can indeed be justified within the 
framework defined by the assumptions summarized at the 
beginning of this section. To illustrate this point, let us 
concentrate first on a]] and Oq inside the ellipse resulting 
from the “scalar” fit. Choosing the scale p = 770 MeV ((5i 
is scale-independent), we write 

(5i = 5ct + 5i, (48) 

where the counterterm part (at the scale p) is: 

(5ct = 127r^^(4 + 44 - 44 - 242 )) (49) 

whereas the loop part 6 i is given in terms of O(p^) con¬ 
stants 4) 4,4) 4- We take [29]: 


4 = -0.4 ±0.6, 4 = 4.3 ±0.1, (50) 

and for 4)4 we use the result of the “scalar” fit repre¬ 
sented by the corresponding ellipse in Fig. 3. This leads to 
the estimate 5^ = 0.0004±0.0083. Using the resonance es¬ 
timates of the r’s [9,33,29] which are obtained according 
to the prescription mentioned above and assuming that 
the typical mass scale relevant in the scalar channel is 
Ms = 1 GeV, one obtains 5 ct ~ 0.00025. The authors 
of Ref. [29] assume that this estimate (fi = —1.5,4 = 
3.2,4 = —4.2,4 = —2.5,42 = —0.7) holds within a fac¬ 
tor of 2. We can even relax this detailed estimate in favour 
of a more crude order of magnitude: 

4 = (47r)4r„(/x) ~ ± . (51) 
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which is expected to hold for 0 (p®) constants describing 
symmetry breaking effects in the scalar channel. The cru¬ 
cial point here is the dependence on the effective scalar 
mass Ms- As long as Ms ~ 1 GeV, 6 ct will remain small: 
adding individual f„ contributions randomly (f„ ~ ± 1 - 8 ), 
one gets in this case dcT ~ ±0.0012. The whole two- 
loop correction 6a is then obtained by adding Sm given 
by Eq. (45) with the estimates of dg and Sct- One ob¬ 
tains da = 0.034 ± 0.010 in agreement with the value of 6 a 
required by the scalar ht and shown in Table 2. 

We now return to the result of our paper indicating a 
deviation of the CGL narrow strip from experiment. Let 
us assume for the moment that the value of the scalar 
radius remains unchanged. Then, according to Table 2, 
the actual value of the two-loop correction 6 a = 0.020 ± 
0.015 is no longer saturated by 6 m = 0.037 ± 0.005 and 
consequently (5i can no longer be negligible. The loop part 
6e can be estimated as before. Taking the same values 
Eq. (50) of ii ,£2 and using for ^ 3,^4 the results of the 
global ht (see the corresponding ellipse on Fig. 3), one 
obtains in this case^ 61 = 0.0096 ± 0.0076. This allows the 
extraction of the required value of the 0 {p^) counterterm 
combination 6ct- Parametrizing the possible variation of 
the scalar radius by: 

(r^), =0.61 fm2(l±,5^), (52) 

one obtains the hnal estimate for the combined effect of 
6ct and 6r (our analysis does not allow us to disentangle 
these corrections): 

6ct + 0.05541 • ,5^ = -0.0266 ± 0.009 . (53) 

This means that both 0{p^) contributions - the subtracted 
infrared singular moment 6m and the symmetry breaking 
counterterms - are basically of the same order of magni¬ 
tude. This is still consistent with the crude order of magni¬ 
tude estimate of Eq. (51), provided the effective mass scale 
Ms characteristic of the scalar channel contributions is re¬ 
duced by a factor 2: Ms ~ 500 MeV. This could indeed 
be a rather natural (though rough and qualitative) way 
how to account for the exceptional role of the tttt contin¬ 
uum and of the OZl-rule violation in the scalar channel. 
In this case, the estimate Eq. (51) leads to |(5ct| ~ 0.019, 
in qualitative agreement with Eq. (53), provided 6ct is 
negative. 


8 Conclusion 

Low-energy tttt scattering has long been recognized as the 
golden observable to access the chiral structure of QCD 
vacuum. From one side the Roy equations allow a com¬ 
pletely model-independent experimental determination of 
the two S-wave scattering lengths Oq and Oq. Chiral sym¬ 
metry can then be used to translate this information into 
knowledge of the EEC’s of the chiral Lagrangian. In this 

® All the above estimates of the loop part of Ai can be re¬ 
produced by using its expression in terms of a% and Uq, which 
arises from the matching with the Roy equations’ solutions. 


paper we have followed these two steps using recent data 
on A'e 4 decays published by the E865 Collaboration. Con¬ 
trary to previous analyses, we did not rely on any theoret¬ 
ical assumption about the correlation between Og and Oq, 
but rather supplemented the data with existing data 
in the 1 = 2 channel below 800 MeV. The result 

o); = 0.228 ±0.012, ag =-0.0382 ±0.0038 (54) 

is compared with the theoretical relation between 2 aQ—5ag 
and the scalar radius of the pion [ 10 ] obtained in the stan¬ 
dard two-loop yPT. If the dispersive determination of the 
latter, {r^)s = (0.61 ± 0.04) fm^ is used, one finds a dis¬ 
agreement at the 1-a level. It is possible to turn the argu¬ 
ment around, and interpret this discrepancy as a measure¬ 
ment of 0 (p®) counterterms that contribute to the above 
theoretical relations. The latter come out larger than ex¬ 
pected by the usual resonance saturation assumptions; 
this fact might be the manifestation of the exceptional 
status of the scalar channel, characterized by a strong tttt 
continuum and OZI rule violation. 

The same conclusions are reached once we include in 
the Roy equations solutions the dependence on the phases 
at the matching point, -^/s = 800 MeV, Oq and 0i. These 
two quantities are likely to represent the most important 
source of theoretical uncertainty for the Roy solutions, 
aside from the one on the driving terms. Our extended 
Roy solutions concretely parametrize the theoretical error 
on the phaseshifts, for a given value of Oq and Oq. 

Comparing the full two-loop standard yPT prediction 
for the two scalar scattering lengths [29], ajj = 0.220 ± 
0.005 and Oq = —0.0444 ± 0.0010 with Eq. (54), one Ends 
agreement for ajj, but a 1 -a discrepancy persists for Oq. 
This makes still more crucial the outcome of new precise 
experiments on tttt, which are either planned [25] or on¬ 
going [24]. In particular, more accurate Ke 4 data in the 
region of higher could eventually allow a simultane¬ 
ous determination of both Og and Og from a single set of 
data on — This might be preferable to the procedure 
adopted in the present paper, where E865 results are com¬ 
bined with the TT+TT^-production data of a different ori¬ 
gin. Our result on ajj Eq. (54) agrees with one of the two 
determinations of the isoscalar scattering length by the 
E865 collaboration [1] that does not use the narrow strip 
correlation between ajj and Og as a theoretical input [ 10 ]. 
No determination of Og is reported in [1]. In this respect 
it is interesting to notice that the observed experimental 
correlation between ajj and Og is positive and close to 1. 
The pionium lifetime experiments [24] cannot distinguish 
(ag,ag) pairs with the same [og — Og], in particular the 
global and scalar Hts which lead to Og-Og = 0.266±0.010 
and 0.263 ± 0.010 respectively. 

After determining the scattering lengths, we have stud¬ 
ied the consequences for the tttt subthreshold parameters, 
by means of a matching procedure with the chiral ampli¬ 
tude. The influence of the uncertainty in 0g and 0i is only 
apparent in the parameters Ai and A 2 , whereas the others, 
in particular a and /3, are practically independent thereof. 
These last two parameters are intimately related to the 
two main order parameters of SU (2) x SU (2) chiral sym¬ 
metry, the quark condensate and the pion decay constant 
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in the SU{2) chiral limit, keeping the strange quark mass 
at its physical value. The result for the two-flavour GOR 
ratio ^(2) = 0.81 ± 0.09 corresponds to a large and nega¬ 
tive central value for the LEG £3 ~ —18, to be compared 
with the standard expectation £3 = 2.9 ± 2.4. Since this 
constant is rather sensitive to the OZI rule violating con¬ 
stants L 4 and Lq of the three-flavour chiral Lagrangian, 
this conclusion seems to corroborate previous indications 
concerning the size of OZI rule violation [12,15] in the 
scalar channel, namely the value of Lq. 

Finally, we would like to stress that, although crucial 
for understanding the pattern of Nf = 2 chiral symmetry 
breaking, tttt scattering alone does not tell us anything 
about the important question of the dependence of chi¬ 
ral order parameters on Nf. In particular it cannot be 
used to disentangle the “genuine” condensate of the purely 
massless theory, from the one “induced” by massive (but 
light) strange quark pairs. Nevertheless very accurate low- 
energy data on tttt scattering still represent an essential in¬ 
gredient in a combined SU{3) x SU{3) analysis and forth¬ 
coming determination [23] of the three flavour condensate 
27(3), the quark mass ratios and other “strange” features 
of the QCD vacuum. For this program, new low-energy 
TTTT and ttK scattering data are of obvious interest. 
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sallam and H. Leutwyler for useful discussions and correspon¬ 
dence, J. Gasser for useful comments on the manuscript. Work 
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0169 (EURODAPHNE). SD acknowledges partial support by 
PPARC, through grant PPA/G/S/ 1998/00530. LG acknowl¬ 
edges partial support from European Program HPRN-CT-2000- 
00149. 


A Treatment of errors for ACM(A) data 

As discussed at the end of Section 3, it was proposed in 
Ref. [5] to take into account (unknown) systematic errors 
in method A of Ref. [ 6 ], by enlarging errors of AGM(A) 
solution in the following way. Consider for each energy bin 
the phase shift of solution A { 6 a ± 2^6a) and the one of 
solution B (Sb^^Sb)- The difference [( 5^1 —fe] is supposed 
to reflect systematic errors of method A. The prescription 
advocated in Ref. [5] consists then in adding in quadrature 
this difference to the error A6 a quoted in Ref. [ 6 ]. The 
ACM (A) solution with enlarged errors is thus defined as: 

6 a ± \/ {A6Ay + {6 a — 6 b)‘^ (55) 

The “global” and “extended” fits can be performed 
with these enlarged errors for ACM(A) phase shifts. The 
corresponding 1 - and 2-a contours for the S-wave scatter¬ 
ing lengths are plotted in Fig. 6 , while Table 3 summarizes 
the changes occurring to the various quantities considered 
in this paper. Let us mention that the nonlinear Eqs. (28), 
(32) and (33) have been used to compute £3 and £4, and 
a contribution of the NNLO remainders 6 ,e < 1% is in¬ 
cluded in the error bars of X{2) and /F^. 


Table 3. Results of the global and extended fits, considering 
for ACM (A) phase shifts either the errors indicated in Ref. [6], 
or enlarged errors according to the prescription of Ref. [5]. 





Global ht 




1 ACM(A) errors 

1 Enlarged errors | 

0.0 

0.228 

± 

0.012 

0.227 

± 

0.012 

al 

-0.0382 

± 

16.45 

0.0038 

-0.0392 

± 

14.56 

0.0041 

a 

1.381 

± 

0.242 

1.334 

± 

0.252 

p 

Pa.13 

1.081 

± 

-0.14 

0.023 

1.088 

± 

-0.25 

0.024 

£4 

-17.8 

± 

15.3 

-14.1 

± 

15.0 

£4 

4.1 

± 

0.9 

4.3 

± 

0.9 

X{2) 

0.81 

± 

0.07 

0.82 

± 

0.07 

{F/F.f 

0.89 

± 

0.02 

0.89 

± 

0.03 





Extended ht 




1 ACM(A) errors 

1 Enlarged errors | 


0.228 

± 

0.013 

0.227 

± 

0.013 

al 

-0.0380 

± 

16.48 

0.0044 

-0.0389 

± 

14.59 

0.0047 

a 

1.384 

± 

0.267 

1.340 

± 

0.281 

p 

Paf3 

1.077 

± 

-0.23 

0.025 

1.084 

± 

-0.30 

0.027 

£3 

-18.5 

± 

16.7 

-15.0 

± 

16.3 

£4 

4.0 

± 

0.9 

4.2 

± 

1.0 

X{2) 

0.81 

± 

0.08 

0.82 

± 

0.08 

{F/F.f 

0.90 

± 

0.03 

0.89 

± 

0.03 


We see that enlarging errors for ACM(A) phase shifts 
slightly moves the global/extended contours in the (oq, Uq) 
plane towards the bottom of the Universal Band. The ef¬ 
fect is however marginal: the l-u discrepancy remains be¬ 
tween the scalar and the global/extended fits, and the 
derived quantities £ 3 , £4, 2 ^( 2 ) and F'^/F'^ are almost un¬ 
changed. 

B Parametrization of the Roy solutions 

In Ref. [5], the solutions of Roy equations have been de¬ 
scribed according to the Schenk parametrization Eq. (5). 
The Schenk parameters Xj {X = A, B, C, D, s) are func¬ 
tions of the scattering lengths, Oq and ag, and the phase 
shifts at the matching point, 6 q{so) = 6*0 and 1/1 (■^o) = 

In Ref. [5], the Roy equations were solved for the particu¬ 
lar choice of phase shifts at the matching point: 6*0 = 82.0° 
and 01 = 108.9°, and the dependence on ajj and Og of the 
parameters X = A, B, C, D, s was parametrized according 
to Eq. ( 6 ). 

We have followed the same procedure as in Ref. [5], 
with the only difference that do and 9i are explicitly treated 
as variables. After generating Roy solutions for dg G {78.9°, 
82.3°, 85.7°} and di G {106.9°, 108.9°, 110.9°}, we then 
performed a fit of the form of Eqs. (5), ( 6 ) and ( 8 ) with 
our solutions, in order to obtain the parameters ai,bi,Ci 
of Eq. ( 8 ). The Schenk parameters Sg, s} and Sg are not 
parametrized in the form of Eq. ( 6 ), but are fixed by the 
condition at the matching point 6 g{so) = 9i. 
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The parameters a* are obtained by considering only 
the solution (^o = 82.3°, = 108.9°). The parameters 

bi have then been obtained by fitting Roy solutions with 
69o 7 ^ 0 and S9i = 0 , and the parameters Ci with 69o = 0 
and 69i ^ 0. At each step, the fit has been performed 
at the level of the phase shifts, and not of the Schenk 
parameters, in order to ensure a smooth dependence on 
(og, Og, 6 * 0 , ^i). We have checked that this parametrization 
was adequate, even for solutions with both non-vanishing 
(500 and 59i. The maximal gap between any solution and 
our parametrization is at most 1 % in the / = 0 , 1 , 2 chan¬ 
nels. 

The coefficients resulting from the fit are expressed in 
units of Mtt) and the phase shifts 0 o,i ,2 are in radians. For 
Aq and Ag, all the coefficients a^, bi, Ci vanish, apart from: 

Aj] : ai = 02 = 0.225, (56) 

Ag : ai = 03 = -0.03706, (57) 


Par . 

-Zi 

Oi 

h 

Cl 


1 

.3617 • 

10 -^ 

-.1716 • 

10-2 

-.3860 • 

10-2 


2 

.1574 • 

10-1 

-.2448 • 

10-2 

-. 3384 • 

10-3 


3 

.1057 • 

10-1 

-.1774 • 

10-2 

-.2510 • 

10-4 


4 

-.1782 • 

10-2 

-.1025 • 

10-1 

-.4312 • 

10-2 


5 

.2572 • 

10-3 

-.4649 • 

10-2 

-.1705 • 

10-2 


6 

-.2872 • 

10-3 

.1046 • 

10-2 

-.3467 • 

10-2 


7 

.8311 • 

10-2 

-.9152 • 

10-2 

-.3637 • 

10-2 


8 

-.2603 • 

10-2 

-.1489 • 

10-1 

.2188 • 

10-2 


9 

.1247 • 

10-2 

.7639 • 

10-3 

-.1340 • 

10-2 


10 

-.1186 • 

10-3 

.4371 • 

10-2 

.1128 • 

10-4 


Par. 

-Zi 

at 

bi 

Cl 


1 

.2482 

.4902 • 10-1 

.1282 • 10-1 


2 

.1997 

.1630 

-.3179 • 10-3 


3 

.1285 

.1137 

.1640 • 10-3 


4 

.1831 • 10-1 

-.1185 

.6305 • 10-1 

rO 

-^0 

5 

.9970 • 10-2 

-.6395 • 10-2 

.1104 • 10-1 


6 

.4846 • 10-1 

.3431 

-.1661 • 10-1 


7 

-.3888 • 10-2 

-.1598 

.4322 • 10-1 


8 

-.8912 • 10-2 

.5183 

-.3067 • 10-1 


9 

-.4265 • 10-2 

.4161 • 10-1 

.8623 • 10-2 


10 

-.3232 • 10-2 

-.1073 

.2976 • 10-2 


Par. 

Zi 

ai 

bi 

Cl 


1 

.1135 • 

10-3 

-.1685 • 

10-3 

-.6043 • 

10-3 


2 

-.2094 • 

10-2 

-.3429 • 

10-3 

-.5583 • 

10-4 


3 

-.8626 • 

10-3 

-.2467 • 

10-3 

-.2205 • 

10-4 


4 

.2911 • 

10-3 

-.8897 • 

10-3 

-.5793 • 

10-3 


5 

.7343 • 

10-4 

-.4099 • 

10-3 

-.2258 • 

10-3 


6 

.2063 • 

10-3 

-.4832 • 

10-3 

-.6376 • 

10-3 


7 

.5294 • 

10-3 

-.6346 • 

10-3 

-.3879 • 

10-3 


8 

-.3372 • 

10-3 

-.2347 • 

10-2 

.9292 • 

10-3 


9 

-.1564 • 

10-3 

.1032 • 

10-4 

-.1169 • 

10-4 


10 

-.1301 • 

10-4 

.8137 • 

10-3 

-.1051 • 

10-3 


Par. 

Zi 

ai 

bl 

Ci 


1 

-.8567 • 

10-3 

-.5496 • 

10-2 

.1526 • 

10-2 


2 

-.1561 • 

10-3 

.1510 • 

10-2 

-.6254 • 

10-3 


3 

-.8722 • 

10-2 

.9679 • 

10-3 

.2538 • 

10-3 


4 

.9872 • 

10-2 

.1001 • 

10-3 

.2140 • 

10-3 

Bl 

5 

.2176 • 

10-3 

.3724 • 

10-2 

.3595 • 

10-2 

6 

.3338 • 

10-3 

-.1050 • 

10-3 

-.5945 • 

10-2 


7 

-.2051 • 

10-3 

.4012 • 

10-3 

.1157 • 

10-3 


8 

-.5171 • 

10-3 

.7078 • 

10-2 

.1593 • 

10-2 


9 

-.5929 • 

10-3 

-.6046 • 

10-2 

.1382 • 

10-2 


10 

-.2247 • 

10-3 

.4017 • 

10-2 

-.1490 • 

10-2 


Par. 

Zi 

tti 

bl 

Ci 


1 

-.1652 • 

10-3 

.2246 • 

10-3 

.3320 • 

10-2 


2 

.3280 • 

10-2 

.5387 • 

10-3 

.9391 • 

10-4 


3 

.1127 • 

10-3 

.2911 • 

10-3 

.2303 • 

10-3 


4 

.1367 • 

10-3 

.1198 


.9361 • 

10-2 

CgO 

5 

.1606 • 

10-3 

.5107 • 

10-3 

.1440 • 

10-3 

6 

.2990 • 

10-3 

-.1170 • 

10-3 

.1345 • 

10-2 


7 

-.5982 • 

10-2 

.9021 • 

10-3 

.1428 • 

10-3 


8 

.1923 • 

10-2 

.9601 • 

10-3 

-.4036 • 

10-2 


9 

.1106 • 

10-3 

.2148 • 

10-3 

-.1501 • 

10-2 


10 

.3809 • 

10-2 

-.2854 • 

10-3 

.2780 • 

10-2 



Bi 

Qi 

bl 

Ci 


1 

-.7257 • 

10-4 

-.1076 • 

10-4 

-.8750 • 

10-4 


2 

.2234 • 

10-3 

-.4577 • 

10-4 

-.8053 • 

10-3 


3 

.3718 • 

10-4 

-.3531 • 

10-4 

-.6497 • 

10-3 


4 

.2259 • 

10-4 

.2031 • 

10-4 

-.7306 • 

10-4 

C\ 

5 

.1216 • 

10-4 

-.2042 • 

10-4 

-.2856 • 

10-4 


6 

.4075 • 

10-4 

-.1625 • 

10-3 

-.1121 • 

10-3 


7 

-.1238 • 

10-4 

-.3676 • 

10-4 

-.2568 • 

10-4 


8 

.1103 • 

10-3 

-.3679 • 

10-3 

-.5010 • 

10-4 


9 

.3813 • 

10-4 

-.5706 • 

10-3 

.3202 • 

10-4 


10 

.3531 • 

10-4 

.1373 • 

10-3 

-.3439 • 

10-4 



Bi 

ai 

bl 

Ci 


1 

-.7557 • 

10 --^ 

.2648 • 

10-2 

-.5166 • 

10-3 


2 

.3425 • 

10-1 

-.2038 • 

10-2 

.5412 • 

10-3 


3 

.2830 • 

10-1 

-.9686 • 

10-3 

.2995 • 

10-3 


4 

.3342 • 

10-2 

.5536 • 

10-2 

-.5538 • 

10-2 

^2 

5 

.1391 • 

10-1 

.7956 • 

10-3 

-.2012 • 

10-2 

6 

.2375 • 

10-1 

.1775 • 

10-2 

.2675 • 

10-2 


7 

-.3024 • 

10-1 

-.1924 • 

10-1 

.6680 • 

10-4 


8 

-.9323 • 

10-1 

-.2108 • 

10-2 

.2173 • 

10-2 


9 

-.8813 • 

10-1 

.4251 • 

10-2 

-.2462 • 

10-2 


10 

-.2679 • 

10-1 

-.3504 • 

10-2 

.1984 • 

10-2 



Bi 

ai 

bl 

Ci 


1 

-.6396 • 

10-3 

.7978 • 

10-3 

.6667 • 

10-3 


2 

-.4143 • 

10-2 

.5649 • 

10-2 

-.5508 • 

10-4 


3 

-.3708 • 

10-2 

.5227 • 

10-2 

.1462 • 

10-3 


4 

-.4016 • 

10-2 

-.6414 • 

10-2 

-.8673 • 

10-3 

^ 0 ° 

5 

-.3159 • 

10-2 

-.3022 • 

10-2 

-.9427 • 

10-3 

6 

-.7352 • 

10-2 

.1584 • 

10-4 

.2274 • 

10-2 


7 

-.1305 • 

10-2 

-.1363 • 

10-4 

.3488 • 

10-2 


8 

-.4523 • 

10-2 

.1960 • 

10-4 

.1146 • 

10-2 


9 

-.4581 • 

10-2 

-.2917 • 

10-3 

-.1778 • 

10-2 


10 

-.1272 • 

10-2 

-.4082 • 

10-2 

.1184 • 

10-2 
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Par. 

Zi 

ai 

h 

Cl 


1 

.6607 • 

Tb^ 

-.1767 • 

10-5 

-.1271 • 

10-4 


2 

-.1750 • 

10 -* 

-.5895 • 

10-5 

-.8847 • 

10-5 


3 

-.6507 • 

10-5 

-.5144 • 

10-5 

-.1517 • 

10-5 


4 

-.3851 • 

10-5 

.1657 • 

10 -* 

-.7559 • 

10-5 

Dl 

5 

.4987 • 

10-5 

.2201 • 

10-5 

-.3089 • 

10-5 


6 

.1953 • 

10-5 

-.3159 • 

10 -* 

-.1827 • 

10-4 


7 

-.2797 • 

10 -* 

.3893 • 

10-5 

.1194 • 

10-5 


8 

.1604 • 

10 -* 

-.5762 • 

10 -* 

-.1570 • 

10-4 


9 

-.1183 • 

10 -* 

-.9919 • 

10-5 

.9930 • 

10-5 


10 

-.7835 • 

10-5 

.2179 • 

10 -* 

-.7949 • 

10-5 


Par. 


ai 

bi 

Cl 


1 

.1980 • 

10-4 

.1510 • 

10-4 

-.2527 • 

10-4 


2 

-.2572 • 

10-2 

-.5907 • 

10-4 

.1149 • 

10-4 


3 

-.2024 • 

10-2 

-.2137 • 

10-4 

.1067 • 

10-4 


4 

.1600 • 

10-2 

.5689 • 

10-4 

-.2189 • 

10-4 

n2 

-^0 

5 

.1790 • 

10-4 

.1280 • 

10-4 

-.7452 • 

10-4 

6 

.1228 • 

10-2 

-.3551 • 

10-4 

.9342 • 

10-4 


7 

.9168 • 

10-4 

-.7961 • 

10-4 

-.1405 • 

10-4 


8 

.4960 • 

10-2 

-.1981 • 

10-4 

.7174 • 

10-4 


9 

.5225 • 

10-2 

.2966 • 

10-4 

-.8173 • 

10-4 


10 

.1550 • 

10-2 

-.1694 • 

10-4 

.6938 • 

10-4 


We determine Sq, sJ and Sq by the conditions: 

'^ o ( so ) = 6 ' o , ' 5 i ( so ) = 6 ' i , ( 5 o ( so ) = 6»2 ( 58 ) 

where ^2(00: ®o> ^O) ^i) is parametrized following Eqs. ( 6 ) 
and ( 8 ), with the coefficients: 


Par. 

2 * 

ai 

bi 

Ci 


1 

-.3160 

.7038 • 10-* 

-.2480 • 10-* 


2 

-.2355 

.2380 • 10-* 

.6701 • 10-2 


3 

-.2021 

.1687 • 10-* 

.5869 • 10-2 


4 

.4885 • 10-* 

.6057 • 10-* 

-.2094 • 10-* 

02 

5 

-.1106 • 10 -* 

.2317 • 10-* 

-.1128 • 10 -* 


6 

.8406 • 10-2 

.7702 • 10-* 

-.2254 • 10-* 


7 

.3569 • 10-2 

.1531 

.1103 


8 

.3021 • 10-* 

.1027 • 10-2 

-.4945 • 10-2 


9 

.2762 • 10-* 

.2859 • 10-2 

-.1297 • 10-* 


10 

.7229 • 10-2 

.1513 • 10-* 

.1340 • 10-* 


We have introduced the function 92 {aQ,aQ,6o,9i) to 
improve the accuracy on Sq. The dependence of 02 on the 
S-wave scattering lengths and on the phases at the match¬ 
ing point is smoother than the analogous dependence of 
Sg. A direct fit of Sg using the form of Eqs. ( 6 ) and ( 8 ) 
would therefore induce a loss of accuracy, compared to the 
procedure we follow here. 
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